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Abstract. Landau, Pollak, Slepian, and Tracy, Widom discovered that cer- 
tain integral operators with so called Bessel and Airy kernels possess commut- 
ing differential operators and found important applications of this phenomena 
in time-band limiting and random matrix theory. In this paper we announce 
that very large classes of integral operators derived from bispectral algebras 
of rank 1 and 2 (parametrized by lagrangian grassmannians of infinitely large 
size) have this property. The above examples come from special points in these 
grassmannians . 



It was discovered by Landau, Pollak, Slepian |2()l 1151 ITfll 119) and Tracy, Widom 
|21l |22] that certain integral operators associated to the Airy and Bessel special 
functions posses commuting differential operators. They found important applica- 
tions of this to time-band limiting, and to the study of asymptotics of Fredholm 
determinants, relevant to scaling limits of random matrix models. We call this 
phenomena the prolate spheroidal phenomena. 

On the other hand, the problem of bispectrality was posed [7] about 20 years 
ago by one of us (A. G.) and J. J. Duistermaat as a tool to understand this prolate 
spheroidal property of integral operators. The aim was to extend it to larger classes 
and search for possible applications. Despite the dramatic recent developments in 
the areas of random matrices and bispectrality the two problems remained isolated 
except for several common examples, see |?H 1101 ITT] . In addition only few integral 
operators possessing the prolate spheroidal phenomena were found. 

Here we announce that any selfadjoint bispectral algebra of ordinary differential 
operators (see Definition ^ . 31 and Definition ^ . 31 and Section 2 for general definitions) 
of rank 1 and 2 induces an integral operator possessing the prolate spheroidal 
property. The kernel of such an operator is of the form 



where ^(x,z) is the corresponding bispectral (eigen)function and is a contour 
in the complex plane with 1 or 2 end points. It acts on the space L 2 (r!) again for 
a contour with the same property. The main results are stated in Theorem 13.81 and 
Theorem Ol 

The integral operators of Landau, Pollak, Slepian, and Tracy, Widom correspond 
to the cases ^(x, z) = ^/xzJ^ + i/ 2 (ixz) and ^(x, z) = A(x + z) which are known to 
be "basic" bispectral functions, in the sense that the other rank 1 and 2 bispectral 
functions are obtained from them by certain type of Darboux transformations, see 
El H LI HI (Here J„(.) denote the Bessel functions of first kind and A(.) 
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denotes the Airy function.) Although in these cases [23 EH EH EH EH EU the 
commuting differential operator is of order 2, in general it is of arbitrarily large 
order. 

In the rest of the introduction we describe our strategy for proving Theorem 
13.81 and Theorem 14.81 which relies on a very interesting property of the "size" of 
bispectral algebras of rank 1 and 2. 

Consider, more generally, a holomorphic function ty(x, z) in some domain of 
C x C which is not an eigenfunction of any differential operator in x or z. Denote 
by By the algebra of differential operators R(x,d x ) with rational coefficients for 
which there exists an operator S(z,d z ) with rational coefficients such that 

(1.1) R(x, d x )V{x, z) = S(z, d a )*(x, z). 

The algebra of all differential operators S(z, d z ) obtained in this way will be denoted 
by Cy . The equality 

by(R(x,d x )) :=S(z,d z ) 

correctly defines an antiisomorphism from By to Cy. Recall that such a function 
ty(x, z) is called bispectral if both algebras By and Cy contain rational functions. 
The subalgebra of By and Cy, consisting of differential operators in x and z for 
which ^(x, z) is an eigenfunction, are commutative. Algebras obtained in this way 
are called bispectral algebras. Their rank (which is equal and is also called rank of 
the bispectral function *ff(x, z)) is the greatest common divisor of the orders of the 
operators of these algebras. 

We derive our main result from the following remarkable property: 

Consider the Z + x Z + filtration of the algebra By given by 

B% M = {R{x,d x ) e By | 0TdR(x,d x ) < 2l u ord(byR)(x,d x ) < 2l 2 }. 

If *&(x, z) is a bispectral function of rank r = 1 or 2 then the the size of the 
spaces B 1 ^ 1 " 2 in this filtration is large in the sense 

dim Si 1 '' 2 > -(2l x l 2 + h + l 2 ) - conts 
r 

where the constant is independent of l\ and l 2 . 

For the basic bispectral functions e xz of rank 1, and A{x + z), ^/xzJ^ + i/ 2 (ixz), 
v E C\Z of rank 2, the dimension of these spaces is exactly equal to the right 
hand side with const = — 1. This is remarkable since it shows that for all bispectral 
functions ^(x, z) of rank 1 and 2 the spaces of the above natural filtration of the 
algebra By are approximately of the same dimension as the basic (Bessel and Airy) 
ones. 

It is very interesting to understand the relation between our results and the 
approach of Adler, Shiota, and van Moerbeke E] to the Tracy- Widom system 
of differential equations via Virasoro constraints. One possible way would be to 
incorporate the representation theoretic meaning of the Bessel tau functions from 
[S] in terms of representations of the Wi +00 algebra, see |S]. 

Another important problem is to understand the relation of this work to the 
isomonodromic deformations approach to random matrices from the works of Palmer 
P| . Harnad-Tracy-Widom jT3H2S], Its-Harnad jE], and Borodin-Deift 0. 

The proofs of the results announced here will appear in a forthcoming publica- 
tion. 
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2. BISPECTRALITY AND COMMUTATIVITY 

Denote by W ra t the algebra of differential operators in one variable with rational 
coefficients. Denote by a(.) the formal adjoint of operators in W ra t '■ 

(n \ n 

Y,h(x)d2 ) =^2(-d x rb k (x). 
k=0 I k=0 

For any £ € C consider the map from the subalgebra of W ra t consisting of operators 
regular at £ to the set of bidifferential operators at £ 



\fe=0 / k=0 i=0 



(2.2) =X)E(- 1 ) < ^" < " 1 ®^ 6 *( a ) x=f 



For an oriented contour T in C we will denote its endpoints by e(r). For any 
£ G e(r) set 7r(£) = 1 or depending on whether £ is a left or right endpoint of T. 
Assume that 



D(x,d x ) = J2b k (x)8ZeWr 



Vrat 

k=0 



is a differential operator which is regular along T. Let f(x) and g{x) be smooth 
functions on V that decrease rapidly when x — ► oo (i.e. ]im x ^ >00 p(x)f( k '(x) = for 
any polynomial p(x) and any integer k, similarly for g(x)). 
By standard integration by parts one gets that 

(2.3) f (D(x,d x )f(x))g(x)dx= £ {-l)*®<l>t(D(x,d x ))(f(x)®g(,x)) 
r «ee(r) 

+ J f(x) (aD(x,d x )g(x))dx. 

Let ^{x, z) be a holomorphic function on some domain in C x C which is not an 
eigenfunction of a differential operator in x or z, as in the introduction. Recall the 
definition i|l.lfl of the algebras of differential operators By and Cy. Assume that By 
and Cy are two subalgebras of By and Cy, respectively, that are stable under the 
formal adjoint map (|2.1|l and such that 

by (By) = Cy. 

Proposition 2.1. Let T± and T2 be two contours in C such that T\ x T2 is in the 

domain of ^f(x, z) and ^(x, z) decreases rapidly when x or z go to 00 along T± and 
T2. If D(x,d x ) € By is such that 

(2.4) aby(D(x, d x )) = bya(D(x, d x )) 
and 



(2.5) ^(D(x j a a )) = 0,Ve€e(T 1 ); &((&!>)(*, d z )) = 0, V£ G e(F 2 ) 
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then the integral operator with kernel 

(2.6) K(x,y)= [ *{x,z)V{y,z)dz 

on L 2 (Ti) commutes with the differential operator D(x, d x ) with domain all smooth 
functions on T\ that decrease rapidly as x — ► oo. 

Let us call a differential operator D(x,d x ) € W ra t formally symmetric if 

(aD)(x,d x ) = D(x,d x ) 

and formally skewsymmetric if 

(aD)(x,d x ) = -D(x,d x ). 

Lemma 2.2. (i) A differential operator D(x,d x ) G W ra t is formally symmetric if 
and only if it has the form 

n 

(2.7) D(x,d x ) = J2 d * c *W d * 

i=0 

for some integer n and some rational functions Ci(x). 
(ii) The operator D(x,d x ) given by Ij2.7|l satisfies 

fc(D(x,d x ))=0 

for some fixed (eCiJ and only if 

(fijc fc )(O = 0ybr*=l,...,n,i = 0,... ) fc-l. 

For a given function *S>(x, z) as before let 2?* jSym be the subalgebra of Bq, con- 
sisting of differential operators R(x,d x ) for which both R(x,d x ) and (6^i2)(z, d z ) 
are formally symmetric. Set also Cq,. sym := 6* . sym )- Define the vector spaces 

(2.8) 

Bv'X* = { R & d *) e B * sym I oid R{x , d x ) < 2/j and ord(6*i?)(z, d z ) < 2l 2 }, 
(2.9) 

4\tm = G C*, SJ , m | OYd(b^S)(x,d x ) < 2h andordS(z,d z ) < 2l 2 }. 

Clearly C l ^ ym = 

From Proposition 12 . II and Lemma \2. 21 we obtain: 

Theorem 2.3. Assume the conditions from Provosition \2.1\ for the the function 
^(x,z) and the contours T± and r 2 . If either of the following two conditions is 
satisfied then the integral operator with kernel l|2.6|l possesses a formally commuting 
symmetric differential operator of order less than or equal to 2l\ and domain - the 
space of rapidly decreasing smooth functions on Y\. 
Condition (i): 

B*!s Vm > Hh + l)e(ri)/2 + h{l 2 + l)e(r 2 )/2, 

Condition (ii): — e(I\) = e(Fi), — e(T2) — ^,{^2), all operators in Bq, iSym are 
invariant under the transformation x 1— > —x, and 

B*tlm > h(h + l)e(ri)/4 + l 2 (l 2 + l)e(T a )/4 
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3. Integral operators associated to selfadjoint Darboux 
transformations of alry functions 

3.1. The Airy bispectral function. Denote by A(x) the Airy function and set 

(3.1) #a(x,z) = A(x + z). 

Recall that A(x) decreases rapidly when x — ► oo in the sector — 7r/3 < argx < 7r/3. 
If L A (x,d x ) denotes the Airy differential operator 

L A (x,d x ) = d 2 x - x 

then ty A (x,z) satisfies 

(3.2) L A {x, d x )V A (x, z) = z^ A {x, z), 

(3.3) d x ^ A {x,z)=d z ^ A {x,z), 

(3.4) xV A (x,z) = L A (z,d z )V A (x,z). 

For shortness denote the algebras By A and Cy A of differential operators with ra- 
tional coefficients associated to the Airy function *£ A (x, z), recall (|l.l(l . by B A and 
C A . It is straightforward to deduce: 

Lemma 3.1. The algebras B A and C A coincide with the Weyl algebra W po i y of 
differential operators in one variable with polynomial coefficients. Moreover the an- 
tiisomorphism b A associated to the Airy function ^> A (x,z), recall l|3.1|) . is uniguely 
defined from the relations 

b A (x) = (L A (z,d z )), b A (d x )=d z , b A (L A (x, d x )) = z. 

3.2. Selfadjoint Darboux transformations from the Airy function. Note 
that 

C[x] =B A D C(x) 

and 

C[L A (x,d x )] =b A 1 (C A nC(z)). 

The set of rational Darboux transformations T> A from the Airy function was de- 
fined in 4^ as the set of functions ^>(x, z) for which there exist differential operators 



(3.5) P(x, d x ), Q(x, d x ) e (Sa)(c[x]\o) = W ra t 
such that 

(3.6) f(L A (x, d x )) = Q(x, d x )P{x, d x ), 

(3.7) (a;, z) = -^P(x, d x )V A (x, z), 



for some polynomials f(t) and p(z). (The polynomial p(z) is included for normal- 
ization purposes only.) The quotient ring of B A by C[a;]\{0} in (|3.5I) is well defined 
since C[a;]\{0} satisfies the Ore condition, see |17| . 

It was also shown in [4lll4| and more conceptually proved in |3] that: 

Theorem 3.2. All rational Darboux transformations from the Airy function ^(x, z) 
are bispectral functions of rank 2. 
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Definition 3.3. Define the set ST) a of sclfadjoint Darboux transformations from 
the Airy function ^(x, z) to consist of those functions ^{x, z) for which there exists 
a differential operator P(x, d x ) € W ra t such that 

(3.8) g(L A (x, d x )f = (aP)(x, d x )P(x, d x ), 

(3.9) *{x, z) = -^P{x, d x )V A (x, z), 

for some polynomial g{t). 

In fact ST) a consists exactly of those ^(x,z) £ T>a for which Q(x,d x ) = 
(aP)(x, d x ) in (|3 . €>|) — (|3 . T|> with an appropriate normalization of the polynomial p(z). 
One can show that as a consequence f(t) is the square of some polynomial g{t), 
compare to l|3~£|l - (f3~3|) . 

3.3. Size of the algebra Ba relative to the antiisomorphism b A - Consider 
the Z + x Z + filtration of the algebra Ba associated to the Airy function $! A (x, z) 
defined analogously to i|2.8|) by 

(3.10) B l X M = {R(x,d x ) e B A \ ovdR{x,d x ) < 2h,ovd{b A R){z,d z ) < 2l 2 }. 

Lemma 3.4. The vector space B^' 12 has a basis that consists of the differential 
operators 

{x m (L A (x,d x )) n \ n<h,m< l 2 } U {x m d x (L A (x,d x )) n \n<h,m< l 2 }. 

Note that the formal adjoint antiinvolution a of W ra t preserves the spaces Bl' 2 . 
Since a 2 = id the space B 1 )^ 2 is the direct sum of the eigenspaces of a with eigenval- 
ues ±1. Denote the eigenvalue 1 subspace of B l l^ 2 by B l )^ 2 ym . For the Airy function 
^ a(x, z) one has ah a = bAa and thus: 

A,sym 



B A tm = Wx,d x ) e B A \oTdR(x,d x ) < 2h, OTd(b A R)(z,d z ) < 2l 2 , 



aR(x,d x ) = R(x,d x ), a(b A R(x,d x )) = b A R(x,d x )}. 
Lemma 3.5. The set of operators 

(3.11) {x m (L A (x,d x )) n + (L A (x,d x )) n x m \n<h,m< l 2 } 
is a basis for the space B\ ' S y m ■ In particular 

dim^ m = (/ 1 + l)(/ 2 + l). 

3.4. Size of the algebra £>* iSym for a selfadjoint Darboux transformation 
from the Airy function, relative to the involution by. Fix an arbitrary 
selfadjoint Darboux transformation from the Airy function "J € ST) a given by 
(E2I)-<E]II> for s om e P{x,d x ) G Wrat and g(t) € C[t}. Let 

(3.12) P(x,d x ) = -^R(x,d x ) 

for some R(x,d x ) € W po iy — Ba and v(x) G C[x\. Set 

(3.13) ordR(x, d x ) = p\ and ord(bAR){x, d x ) = p 2 . 
Denote 



5*.! = Sptm{-^-R{x,d x )M(x,d x )(aR)(x,d x )-^- \ M(x, d x ) e B^J 2 }, 

V(X) V(X) n,t>u<H 

5*, 2 = Sp<m{v(x)M(x,d x )v(x) | M(x, d x ) e B 1 ^ 2 }. 
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Proposition 3.6. In the above setting: 

(i) The spaces of differential operators <S^ i and 5^,2 are subspaces ofB^'J 1 . 

(ii) The dimension of the intersection S^^± D 5^,2 is less than or equal to 

(h-Pi + l)(h-p2 + l)- 

Theorem 4.2 of [H] shows that (|3.8|1 - 1|3.9(1 is equivalent to 

v{L A (z,d z )f = {ab A R){z,d z )-^(b A R)(z,d z ), 

$(x,z) = 1 (b A R)(z,d z )y A (x,z). 
v{x)g{z) 

(The hard step is to prove the first equality.) One can show that 

6* (-^-R(x,d x )M(x,d x )(aR)(x,d x )-^-) = g(z)(b A M)(z,d z )g(z) 
\v{x) v(x) J 

for all operators M{x, d x ) £ B A . The fact that 6>* ; 2 is a subspace of B l ^' l s 2 ym follows 
from this and the fact that Sis, A C B^'^ ym by exchanging the roles of x and z. 

Theorem 3.7. For any self adjoint Darboux transformation from the Airy function 
~$?(x,z) € ST> A the dimension of the space of differential operators S^i + is 
greater than or equal to (ii + 1)(^2 + 1) — PiP2- In particular 

(3-14) dim£^ m > (h + l)(h + 1) - P1P2- 

Theorem l3.7l and Theorem l2 . 3l implv our final result for integral operators derived 
from Darboux transformations from the Airy function. 

Theorem 3.8. Let ^(x, z) S ST> A be a selfadjoint Darboux transformation from 
the Airy function, given by (|3.8|l . (|3.9[l . (|3.12() . Let T±, T 2 be two connected contours 
in C that do not contain the roots of the polynomials v(t) and g(t) respectively and 
that begin at some finite points and go to infinity in the sector — 7r/3 < argx < 7r/3. 
Then the integral operator 



K(x,y) = / ^(x,z)^(y,z)dz 
Jr? 



on L 2 (Ti) commutes with a formally symmetric differential operator with rational 
coefficients of order less than or equal to 1p\pi and domain all smooth functions 
on Ti that decrease rapidly as x — > 00. 

4. Integral operators associated to selfadjoint Darboux 
transformations of bessel functions 

4.1. The Bessel bispectral function. Denote by J v the standard Bessel func- 
tions of first kind. By abuse of notation the functions 

(4.1) $v(x,z) = (xz) 1/2 J v+1/2 {ixz) 

will be also called Bessel functions. Consider the Euler operator 

D x — xO x 

and the operators 

(4.2) L v (x, d x ) = d 2 x - Ot+H = ±(D X + v){D x -v-1) 

x z x z 



8 



F. ALBERTO GRUNBAUM AND MILEN YAKIMOV 



to be called Bessel operators. The Bessel functions satisfy the equations 

(4.3) L„(x, d x )V v (x, z) = z 2 ty„(x, z), 

(4.4) D x V v (x, z) = D z V v {x, z), 

(4.5) x 2 V u (x, z) = L v (z, d z )$„(x, z). 

For shortness the algebras B^ v and associated to the Bessel function ^(x, z), 
recall Ql.l|) . will be denoted by B v and C„. 

The Bessel functions corresponding to vi,^ £ C that differ by an integer can 
be obtained by a Darboux transformation from each other: 

z) = —{D x -u- z), 

xz 

®»(x, z) = — (D x + u + l)V v+l (x, z), 

xz 

which corresponds to the factorizations 

L v = x~ l {D x + v + ljaT^Ac - v - 1), L v+X = x(D x -v- lJaT 1 ^ + 1/+ 1). 

According to l|4.3[) - (|4.5[) the algebras of differential operators B v and C v contain 
the operators L v (x,d x ), D x , and x 2 . Denote their subalgebras generated by those 
operators by B v and C u , respectively. Clearly H4.3JI - (|4.5JI imply 

b u {B v )=C v . 
Similarly to Proposition 2.4 in J] one shows: 

Lemma 4.1. (i) For v S C\Z i/ie algebras B v andC v are generated by the operators 
L v (x,d x ), D x , and x 2 , i.e. 

By — Bii ; (Lii — Ci/ . 

(ii) For the subalgebras B v and C v of B v and C v consist of exactly those 

differential operators in B v and C„ that are invariant under the transformation 
x i ► —x. 

Note that 

(4.6) ab„(P(x, d x )) = b v a(P(x, d x )), VP(x, d x ) e B„. 

4.2. Selfadjoint Darboux transformations from the Bessel functions. Sim- 
ilarly to the Airy case we have 

C[x 2 } =B V C\ C(x), 
C[L v {x,d x )]=b- x {C v nC{z)). 

The set of rational Darboux transformations X>„ from the Bessel function ^(x, z) 
is defined to be the set of functions ^(x, z) for which there exist differential opera- 
tors 

(4.7) P[x, d x ), Q(x, d x ) e (B V ) (C[X 2 M0}) 
such that 

(4.8) f{L v (x, d x )) = Q(x, d x )P{x, d x ), 

(4.9) tf (i, z) = -jvPti, d w )V u (x, z), 



THE PROLATE SPHEROIDAL PHENOMENA AND BISPECTRALITY 



9 



for some polynomials f(t) and p(z). Again the subset C[x 2 ]\{0} of B v satisfies the 
Ore condition and the quotient ring in (|4.T|> makes sense. 

The following theorem was proved in [21] for v = 1 and in 0] |3j in the general 
case. 

Theorem 4.2. All rational Darboux transformations from the Bess el functions are 
bispectral of rank 2 if v £ C\Z and of rank 1 if v £TL. 

Definition 4.3. Define the set of selfadjoint ("even, selfadjoint" in the case v £ Z) 
Darboux transformations ST) V from the Bessel functions fy u (x,z) to consist of all 
functions ^(x,z) for which there exists a differential operator 

P(x,d x ) e (C[x 2 ]\{0})- 1 ^ 

such that 

(4.10) f(L v (x,d x )) = (-l) m (aP)(x,d x )P(x,d x ) 

(4.H) 9(x, z) = 1 P(x, d x )* u {x, z), 

z m g\z z ) 

for a polynomial f(t) of the form 

(4.12) f(t) = t 2m g(t 2 ) 2 , 5 (0) + 0; g(t) £ C[t]. 

Here the term "even" reflects the fact that for v £ Z the algebras B v and C v are 
bigger than B v and C„. The reason for this terminology is explained below. 

As in the Airy case the set ST) v consists of those rational Darboux transforma- 
tions *&{x, z) from the Bessel function ^ u (x, z) for which in the notation (|4.8|) ~ (|4.9|) 

(4.13) Q(x,d x ) = (-l)° TdP (aP)(x,d x ) 
with the additional property in the case v£Z 

(4.14) P(-x,-d x )=P(x,d x ). 

As a consequence it is obtained that the polynomial f(t) in l|4.8|l has the form (|4.12|) 
and an appropriate normalization of p(z) is made. 

Note that Lemma 14.11 implies that in the rank 2 case v G C\Z the condition 
(|4.14|l is a consequence of l|4.7|l . In the rank 1 case v G Z the term "even" reflects 
this extra condition. It is needed since in the case of Darboux transformations from 
the Bessel function the prolate spheroidal property will be deduced from the second 
condition in Theorem 12.31 

4.3. Size of the algebras B v relative to the antiisomorphisms &„. Consider 
the Z + x Z + filtration of the algebras B v 

(4.15) B l u lM = {R(x,d x ) e B v | oidR(x,d x ) < 2h,0Td(b u R){z,d z ) < 2l 2 } 

The formal adjoint involution a of W ra t preserves the spaces Bj' 2 . Similarly to 
Lemma 13.41 one shows: 

Lemma 4.4. The vector space Bj' * has a basis that consists of the differential 
operators 

{x 2m {L„{x,d x )) n | n < h,m < 1 2 }U {x 2m D x {L v (x,d x )) n \n<h,m< h}. 
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Since a 2 = id the space B l J' 12 is the direct sum of the eigenspaces of a with 
eigenvalues ±1. The eigenvalue 1 subspace of Bl 1 ' 12 will be denoted by B l J t '^ m . The 

commutativity 14.tif) of a and b u on B v implies: 

Bliim = {R(x,d x ) G B v \ ordR{x,d x ) < 2l 1} ord(b„R)(z,d z ) < 2l 2 , 

aR(x, d x ) = R(x, d x ), a{b v R{x, d x )) = b v R{x, d x )}. 

Lemma 4.5. The set of operators 

(4.16) {x 2m (L v (x, 8 x )) n + (L u (x, 8 x )) n x 2m \n<h,m<l 2 } 
is a basis for the space B\)j 2 m . In particular 

dimB l jj 2 m = (h + l)(l 2 + l). 

4.4. Size of the algebra B^ sym for an (even) selfadjoint Darboux transfor- 
mation from a Bessel function. Fix an arbitrary selfadjoint (and in addition 
even in the rank 1 case v G Z) Darboux transformation from a Bessel function 
V v (x,z), W G SV V given by ED} for some P(x,d x ) G (C[a; 2 ]\{0})- 1 6 ly 

and f{t) = t 2m g{t 2 ) 2 , g(t) G C[t\. Let 

(4.17) P{x,d x ) = -^rR{x,d x ) 

for some R(x,d x ) € B v and v(x) G C[x]. Set 

(4.18) ordi?(x, d x ) = p\ and ord(b ll R)(x, d x ) = p 2 . 
Denote 

5*,! = Span{^i?( a; ,^)M( a; ,^)(ai?)(x,9 x )^ | M(x,d x ) G B?,™ h }, 

5*. 2 = SpanMa^AffoQcMz 2 ) | Affofl,) G B^*}. 

Proposition 4.6. (%) TTie spaces of differential operators S^_\ and S^_ 2 are sub- 
spaces ofBy* ym and are invariant under the transformation x \— > —x. 

(ii) The dimension of the intersection S^^\ PI S^_ 2 is less than or equal to 

(h- Pi + l)(l 2 - p 2 + l). 

Finally we obtain the following Theorem. 

Theorem 4.7. For any selfadjoint (and even, selfadjoint in the case Dar- 
boux transformation from a Bessel function ^ v (x,z), ty(x,z) G ST> V the space 
+ <S*,2 consists of differential operators invariant under the transformation 

x i ► —x 

dim(5*4 + Sy, 2 ) > (h + l)(l 2 + 1) - pip 2 . 

In particular the dimension of the subspace of B 1 ^ 'g ym of differential operators in- 
variant under x <— > —x is greater than or equal to (l\ + 1)(^2 + 1) — P\Pi- 

Theorem l4.7l and Theorem l2 . 3l inrplv our final result for Darboux transformations 
from the Bessel functions: 
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Theorem 4.8. Let ^(x,z) G ST> V be a selfadjoint (and in addition even if v G Z) 
Darboux transformation from the Bessel function '^ u {x 1 z) 1 given by (14.10(1 . 1(4.11(1 . 
(14.12(1 . 1(4.17(1 . Let Ti, T2 be two connected finite contours that do not contain the 
roots of the polynomials v(t) and g(t), respectively and such that e(Ti) = — e(Ti). 
Then the integral operator with kernel 

K(x,y) = / z)^(y, z)dz 

on L 2 (Ti) commutes with a formally symmetric differential operator with rational 
coefficients of order less than or equal to 2pip 2 and domain all smooth functions 
on T\. 

Remark 4.9. Let $i(x,z) and $2(2;, z) be two functions satisfying 1(3.2(1 - ((3.4(1 or 
((4. 3(1 - ((4. 4(1 . E.g. in the rank 1 case (section 4, v = 0) one can choose 

<&i{x,z) = e xz , $ 2 (x,z) = e~ xz 

and in the rank 2, Bessel case (section 4) 

$i(x, z) = ^ v {x, z), $ 2 (x, z) = ^ v {x, -z). 

Fix a differential operator P(x, d x ) and a polynomial g(z) that define an (even) self- 
adjoint Darboux transformation from the Airy or Bessel functions as in Definition 
13.31 and Definition ^. 31 Define 

$i(x, z) = — 5— P(x, d x ), and $i(x, z) = P(x, d x ), i = 1, 2 

g{z) z m g{z 2 ) 

in the Airy and Bessel cases, respectively. Then Theorem l3 . 81 and Thcorcm l4.8l hold 
for the integral operator with kernel 

K(x,y)= / y 1 (x,z)V 2 (y,z)dz 

and exactly the same commuting differential operator that appears in those Theo- 
rems. 
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